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Abstract 



We have studied two-body charmless hadronic decays of B mesons into the 
final states (f)K and <j>K*. Using 9.7 million BB pairs collected with the 
CLEO II detector, we observe the decays B~ — > (pK~ and B° — > cfiK* with 
the following branching fractions: B(B~ — > 4>K~) = (5.5 jlf's ±0.6) x 10~ 6 and 
,8(5° -> = (11.5±|;|±l;f) x 10" 6 . We also see evidence for the decays 

5° -» 4>K° and #~ -> ^X* - . However, since the statistical significance is 
not overwhelming for these modes we determine upper limits of < 12.3 x 10 -6 
and < 22.5 x 10~ 6 (90% C.L.) respectively. 
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The phenomenon of CP violation can be accommodated in the Standard Model (SM) 
by a complex phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark- mixing matrix [0. 
Whether this phase is the only source of CP violation in nature remains an open experimental 
question. Studies of charmless B meson decays will certainly play an important role in 
constraining the CKM matrix and testing the SM picture of CP violation. 

Flavor Changing Neutral Currents (FCNC) are forbidden to first order in the SM. Second 
order loop diagrams, known as penguin and box diagrams, can generate effective FCNC 
which lead to b — > s transitions. These processes are of considerable interest because they 
are sensitive to V ts , the CKM matrix element which will be very difficult to measure in 
direct decays of the top quark. They are also sensitive to non-Standard Model physics [0, 
since charged Higgs bosons, new gauge bosons or supersymmetric particles can contribute 
via additional loop diagrams. 

Among final states produced by the gluonic penguin, b — > sg, decay modes in which a 
gluon splits into two strange quarks, g — > ss, play a special role since they cannot be produced 
by any other b decay with comparable rate, thus providing an unambiguous signature for the 
gluonic penguin. As illustrated in Fig. [TJ, a particularly clean final state is produced when 
the kaon includes the spectator quark and no pions are emitted. In this Letter we present 



the first significant measurements of exclusive charmless hadronic decays B — ► <j)K and 
B — > <f>K*. Measurements of other charmless hadronic decay modes by CLEO are discussed 
elsewhere ||. 

The data set used in this analysis was collected with the CLEO detector at the Cornell 
Electron Storage Ring (CESR). It consists of 9.13 fb _1 taken at the T(4S) (on-resonance), 
corresponding to 9.66M BB pairs, and 4.35 fb _1 taken below BB threshold, used for con- 
tinuum background studies. 

CLEO II is a general purpose solenoidal magnet detector ||. Cylindrical drift cham- 
bers in a 1.5T solenoidal magnetic field measure momentum and specific ionization (dE/dx) 
of charged particles. Photons are detected using a 7800-crystal CsI(Tl) electromagnetic 
calorimeter. In the CLEO II. V detector configuration, the innermost chamber was replaced 
by a three-layer, double-sided silicon vertex detector, and the gas in the main drift cham- 
ber was changed from an argon-ethane to a helium-propane mixture. As a result of these 
modifications, the CLEO II. V portion of the data (2/3 of the total) has improved particle 
identification and momentum resolution. 

In the analysis presented here we search for B meson decays by selecting and K (K*) 
decay candidates using specific criteria. Track quality requirements are imposed on charged 
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FIG. 1. Penguin diagram describing B — > cf)K^*' decays 
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tracks, and pions and kaons are identified by dE/dx. Electrons are rejected based on dE/dx 
and the ratio of the track momentum to the associated shower energy in the Csl calorimeter; 
muons are rejected based on their penetration depth in the instrumented steel flux return. 
The Kg candidates are selected from 7t + tt~ pairs forming well-measured secondary vertices 
with invariant mass within three standard deviations (a) of the nominal Kg mass and a 
decay path significance of at least 3<r. The neutral pion candidates are formed from pairs of 
isolated photon-like energy clusters in the Csl calorimeter with invariant mass within —3.5 
and +3.0 standard deviations of the 7r° mass. The meson candidates have K + K~ mass 
within ±20 MeV/c 2 (±4.5T, V = natural width) of the known <ft mass, and the specific 
ionization of the tracks are consistent with the K + K~ hypothesis. The K* candidates are 
reconstructed in four modes: K*° — > K~tt + , K*° — > K°tt , K*~ — > K~7T° and K*~ —> K°7r~, 
and their masses lie within ±75 MeV/c 2 (±1.5r) of the respective known masses. Charmless 
two-body B decays produce the fastest secondary particles among all B decays. Therefore, 
to reduce combinatoric backgrounds only the fastest <fi and the fastest K (K*) are used in 
those events with multiple combinations. 

The B decay candidate is identified via its invariant mass and its total energy. We 
calculate a beam-constrained B mass Mb = y E^ egxa — p B , where p B is the B candidate 
momentum and -Ebeam is the beam energy. The resolution in Mb is dominated by the beam 
energy spread. We define AE = E 1 + E 2 — E hearn , where Ei and E 2 are the energies of 
the daughters of the B meson candidate. We accept events with M B above 5.2 GeV/c 2 and 
\AE\ < 200 MeV. 

We have studied backgrounds from b — > c decays and other b — > u and b — > s decays and 
find that all are negligible for the analyses presented here. The main background arises from 
e + e~ — > qq (where q = u,d,s,c). Such events typically exhibit a two-jet structure and can 
produce high momentum back-to-back tracks in the fiducial region, while BB events tend 
to have more spherical structure, since the B mesons are produced nearly at rest. To reduce 
contamination from these background events, we require the event to have H 2 /H < 0.6, 
where Hi are Fox- Wolfram moments |J. 

We extract the signal yields from unbinned, extended maximum-likelihood fits of the 
preselected on-resonance data separately for each topology (4>K~ , (pK°, (fiK^ K _ n+ , (pK^ KOiT o, 
^iT^-^o, $K*^k°-k-)- F° r an modes, we distinguish signal from background using Mb, AE, 
| cos 6u | (the angle between the thrust axes of the B candidate and that of the rest of the 
event), | cos^l (the angle between the B candidate momentum and beam axis) and (the 
mass of the <p candidate). In addition, we include | cos6^| (the cf> helicity angle, defined as 
the kaon direction in the </> rest frame) for <pK~ and 4>K° modes, irtx* (the mass of the K* 
candidate) for cf)K* modes, and Sk (the number of standard deviations from the predicted 
dE/dx value) of K~ when applicable. 

In each of these fits, the likelihood of the event is the sum of probabilities for the signal 
and background hypotheses, with relative weights determined by maximizing the likelihood 
function C. The probability of a particular hypothesis is calculated as a product of the 
probability density functions (PDFs) for each of the input variables. The signal PDFs are 
represented by a double Gaussian for Mb and AE, by a Breit-Wigner function on top 
of a linear polynomial for and rriK*, by 1 — |cos#b| 2 for |cos#b|, by a third order 
polynomial for | cos 6^1, and by cos 2 6^ for |cos^|. The background distributions for the 
intermediate resonance masses are parameterized by the sum of a Breit-Wigner and a low- 
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order polynomial. For M B , we use an empirical shape (f(M B ) oc M B \/1 — x 2 exp[— 7(1— x 2 )]; 
x = M B /E he£im ) 0. The background AE and | cos^l PDFs are both linear functions, and 
I cos^t| is parameterized by the sum of two terms | cos6 tt \ a with different exponents. The 
signal and background dependences of Sk are bifurcated Gaussian functions. 

The parameters for the PDFs are determined from off-resonance data (background) and 
from high-statistics Monte-Carlo (MC) samples (signal). In the signal MC data, we model 
B — y (j)K^*> as a two-body decay, where for (pK* we assume equal amplitudes for longitudinal 
and transverse polarizations. Dependence on the unknown decay polarization is included in 
the systematic uncertainty. We use a Geant based simulation to model the detector 
response in detail. Further details about the likelihood fit method can be found in Ref. ||. 

To illustrate the fits, we show in Fig. |2] M B and AE projections for the modes with 
significant signals. Events entering these plots must satisfy a requirement on the signal- 
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FIG. 2. M and AE plots for (a,b) B~ 4>K~, and (c,d) B° 4>K*° after the requirement on 
R as described in the text. The projection of the total likelihood fit (solid curve) and the continuum 
background component (dashed curve) are overlaid. For B° — > cj)K* , the two decay modes of K*° 
were weighted according to the statistical errors of the fits. 



to-background likelihood ratios R = \og(P si / Pbi), where P si (Pbi) are signal (background) 
likelihoods for each event i, computed without M B and AE, respectively. This additional 
cut accepts about 2/3 of the preselected signal events in the MC sample. 

We summarize the results for all B decay modes, corresponding submodes, and the 
combined modes in Table |, where we assume equal branching fractions for charged and 
neutral B meson decays 0. We combine the samples from multiple secondary decay channels 
by adding the — 21og£ functions of the branching fraction. The statistical significance of a 
given signal yield is determined from the change in —2 log£ when refit with the signal yield 
fixed to zero. The largest contributions to the systematic uncertainties come, with about 
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TABLE I. Intermediate fitted branching fractions {But), final branching fraction results (B), 
and theoretical estimates [10] are given in units of 10~ 6 . When the result is not statistically 
significant, the final result is quoted as a 90% C.L. upper limit. The errors on branching fractions 
are statistical and systematic respectively. Reconstruction efficiency £ does not include branching 
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ir + 7r decay. 
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i i 9 +3.6+1.8 
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equal weight, from uncertainties in the parametrization of the PDFs, decay polarization 
dependence,^] and possible background from other B decays. 

We observe a significant signal (above 5a) for the decays B~ — > <f)K~ and B° — > (j)K*°. 
Since the statistical significances for the B° — > <pK° and B~ — > <f)K*~ modes are not large 
(2.9a and 3.1cr respectively), we calculate 90% confidence level (C.L.) upper limits (UL) by 
integrating the likelihood curve to 90% of its total area, and increasing it by one unit of the 
systematic error. 

In summary, we have measured B{B- <pR-) = (5.5±?;| ± 0.6) x 10~ 6 and B{B° -> 
(j)K*°) = (11.5^3j1l\j) x 10~ 6 each with statistical significance above 5a. The statistical 
significance of the B° — > <pK° and B~ — > <j)K*~ signals are 2.9a and 3.1a respectively. 
The measured rates are B(B° -> <pK Q ) = (5.4 ± 0.7) x 10" 6 and B(B~ ->■ 4>K*~) = 
{lO-ft-ts-i e) x 10~ 6 . Since the statistical significance in these modes is limited we set upper 
limits of < 12.3 x 10~ 6 and < 22.5 x 10~ 6 (at 90% C.L.) respectively. Averaging over 
B° and B~ we obtain B(B -> <j>K) = (5.5±i;| ± 0.7) x 10~ 6 (6.1a) and B(B -> <f>K*) = 
(11.212 a-ij) x 10 6 (5.9a). The measured branching fractions lie in the range of theoretical 
predictions (see Table [[]). Since there is a considerable spread in theoretical predictions 



1 The B — ► cj)K* decay may be longitudinally or transversely polarized. Assuming 100% transverse 
polarization we obtain B(B — > 4>K*) = (13. 6+44) x 10~ 6 (statistical errors only) and B(B~ — » 
4>K*~) = (12.8±£jj) x 10 -6 . Assuming 100% longitudinal polarization we measure B{B — > 4>K*) = 
(9.9jl|;2) x 10" 6 and B(B~ -» <f>K*~) = (9.9±f;^) x 10~ 6 . To estimate the uncertainty due to the 
unknown polarization we assumed that any value between 100% longitudinal and 100% transverse 
polarization is equally likely. 
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among various calculations, our results will help constrain model parameters. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, the Natural Sciences and 
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Advanced Research Program, and the Alexander von Humboldt Stiftung. 
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